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Abstract 

The representation of the Hubbard operators in terms of the spin— ^ operators 
and the fermion operator with spin— | is proposed. In the low-energy limit this rep- 
resentation is reduced to the representation following from the Hubbard diagramm 
technique. In framework of this approach motion of a hole in an antiferromagnetic 
state of the t-J model is considered. It is shown that the primary hole energy is 
strongly renormalized and the band width has an order of J rather than t. The 
functional integral for the strongly correlated model induced by the obtained repre- 
sentation is formulated. The representation of the total Hubbard algebra for states 
in the lower and the upper Hubbard bands is formulated in terms of the spin— ^ 
and two fermion fields with spin— ^ is formulated. 
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I. Introduction 



The strong correlation in the electron system can be described by the restriction of the 
number of states at every lattice site. Such restriction can be described by the introduction 
of the Hubbard operators. But the Hubbard operators are not convenient for description 
of the statistical and dynamical properties of the strongly correlated electron system. 
That is why some representation of the Hubbard operators in terms of more convenient 
operators are used. The well known representation of such type is slave-boson (fermion) 
representation |l|, ||, |]. These representations are based on the mapping of three electron 
states |0>,|t>)|i>cLta lattice site into fermion and boson Gilbert space and restriction 
of the one-particle states. 

This work has mainly methodical character. At first the simple historical review of the 
description of spin system on the basis of the Holstein-Primakoff and the Dyson-Maleev 
representation is presented. Further the mapping of three Hubbard states into the eight 
spin— I and fermions (with a spin— |) Gilbert space is produced. Such mapping generates 
the representation of the Hubbard operators in terms of fermion and spin— | operators. 
Such representation can be useful in some applications. 

Renormalization of the hole energy in a antiferromagnetic state of the t — J model 
is considered in the next section. This consideration is based on the approach to an 
antiferromagnetic state developed in work [0]. We will show that the method developed 
by Kane, Lee and Read for the description of the motion of a single hole in quan- 
tum antiferromagnet is also applicable to the our approach. We will show that the hole 
spectrum is strongly renormalized in the limit t J where t and J are the hopping and 
exchange constants for the t — J model. The quasiparticle band has the width J and the 
quasiparticle residue is equal to J/t. 

After that two methods of obtaining the functional integral for strongly correlated 
electron system is discussed. In final part of the work the representation of the total 
Hubbard algebra is constructed and its physical interpretation in terms of the local spins, 
electrons in the upper Hubbard bands and holes in the lower Hubbard bands is presented. 

II. Spin systems 

For the clarification of the structure of representation for the Hubbard model we want 
to remember connection between the Holstein-Primakoff (HP) 0] and the Dyson-Maleev 
(DM) Q representations. They based on mapping of the 2S+1 spin states into the 
lowest 2S+1 states of the Heisenberg algebra: 

S-p = P{a^a)^2S ~a+a a P{a+a) 
S+p = P{a^a)a+V2S - a+a P{a+a) 

SIjp = P{a+a){-S + a+a)P{a+a) (2.1) 
here S^,S^,S^ are spin operators, a^,a are creation and annihilation Bose operators. 
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P(a+a) is the projector operator on the lower 2S+1 states: |0 >, | t>,...|2S' >. The 
representation ( |2.1| ) is the direct operator identity: the matrix elements of the original 
spin operator S and the HP spin operators Shp are the same on the lowest 2S+1 states 
of the Heisenberg algebra. The DM representation for spin operators has the form: 

Sdm = -S + a+a (2.2) 

This representation does not represent a direct operator identity but the spin algebra 
is valid for DM spin operators ( |2.2| ) also. One can check that there exists the operator 
l/(a+a) with property 

Shp = PV-'SdmVP (2.3) 

here Shp, Sdm are the HP ( |2.1| ) and the DM ( pl2| ) spin operators. The identity ( |2.3| ) 
means that matrix elements of both parts of Eq.( |2.3|) between any states < n\ and \n > 
are equal to each other. If \n,n > 2S these matrix elements are equal to zero. The 
explicit form of the matrix elements V is following 

n 

< n\V\n >= 5n\n n + 1 - m, n,n <2S (2.4) 

r?i=0 

If we want to calculate the partition function 

Z = Tr{exp{-f3H{S)) = Tr{P exp{-(3H{SHp)) (2.5) 

we can use the identity (|2.3| ) and the commutation of the operators V and P. As a result, 
we have for Z 

Z = Tr{Pexp{-(3H{SDM)). (2.6) 

Eq.( |2.6|) is valid because the action of Sdm on the state \n > for n = 0, 1, 25' leads only 
to the same states, i.e. the lower 25* + 1 states form an invariant subspace with respect 
to the DM spin operators Sdm- 

If we are interested in the properties of the partition function ( |2.5| ) at the low tem- 
peratures in the ferromagnetic or antiferromagnetic state, we can omit the projector P in 
formula (|2.6| ) because the contribution to the higher states \n > with n > 2S into the trace 
( p.6[ ) is exponentially small over the parameter j3J where J is an exchange integral. This 
discussion explains the correctness of using the DM representation for the description of 
the low-energy processes at the low temperatures in ferromagnets and antiferromagnets. 
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III. Strongly correlated systems. 



A similar interpretation can be given for the representation for the Hubbard operators 
X"-^. Let us consider the Gilbert space at every lattice site representing the direct product 
of a Gilbert space of a spin— | and fermion with a spin— |. The total number of states 
in that spin-hole Gilbert space is equal to eight: |0,(t >, \l,a\a >, \2,a >, here the 
first index represents the number of fermions; a, a =t, [ are spin projections. One can 
introduce the Fermi operators h'^, h^r then 

|l,a',cr >= /i+|0,a >, |0, (7 >= /i^' |1, a', cr > . (3.1) 

The singlet and triplet states can be formed from one-fermion states: 

\s >= {i/V2){\i t,i> -|i i,T>), \ti >= |i T,T> 

\to >= (l/v^)(|l T, i> +|i i, T>), |t - 1 >= |i i, i> . (3.2) 

We can map the initial Gilbert space of the Hubbard model without two-fermion states 
into the spin-hole Gilbert space: 

|0>^|s>, \l,a>^\0,a>, (3.3) 

and obtain the following representation of the Hubbard operators X"^ in terms of the 
holes , ha and the spin— | operators s: 

X'^ifp = 2aV2(/i+5)_,(l - n) 
X^^P = 2aV2{l - n){Sh)^a 
Nhp = Xj]p + Xjlp = l-n + d 
X^°p = (1/4) (n - 2d-2{h+ah)s) 

S = s(l - n + d) (3.4) 

here 

5 = (1/4)(1 - 2s(t), n = {h+h), d = h"{h^h+h^. (3.5) 

This representation is the direct operator identity: (1) the matrix elements X'^p operators 
( p.4|) between the physical states |s >, |Ocr > are the same as for the initial Hubbard 
operators X"^; (2) the matrix elements between the unphysical states \tm >, |2cr > are 
equal to zero; (3) the matrix elements between physical states |s >, |0(T > and unphysical 
states \tm >, \2a > are equal to zero. The representation (|3.4| ) for X^fp certainly is 
hermitian and does not require any constraints. If we compare the representation ( |3.4| ) 
with initial expression of the Hubbard operators X"-'^ in terms of the physical Fermi 
operators of electrons we can see that the number of electrons is close to unity while 
the number of Fermi holes is small at filling closed to unity. The empty space without 
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electrons can be imaged as the bound singlet state of the Fermi hole and the spin— | with 
the total spin equal to zero. 

In addition to representation ( |3.4|) in work was proposed the nonhermitian repre- 
sentation for the Hubbard operators 

X^Ta, = 2a/i+(l/2 - (as) + X^^., = 2ah_„, 

NDM = X^j}M + X\,\t = l-n, XOV = n, S = s + (l/2)(/i+a/i). (3.6) 

This representation does not present direct operator identity but the Hubbard algebra is 
valid for this representation. The representation ( |3.6|) can be named as the Dyson-Maleev 
representation for the Hubbard operators. 



What is the connection between representations ( |3.4|) and (|3.6| ) for the Hubbard op- 
erators ? One can check that if we produce canonical transformation of the hole Fermi 
operators in (^]6|) \Pih„, (l/v^)/i+ then the matrix elements of the Hub- 



bard operators in the representation ( |3.4D and ( |3.6| ) are the same between the physical 
states |s >, |Ocr >. Moreover the action of the Hubbard operators in the form (|3.6|) on 
the physical states |s >, \0a > does not lead to the unphysical states \tm >, \2a >, i.e., 
the physical states form invariant subspace over the algebra ( |3.6| ) . The relation between 
the Hubbard operators X^^p (gj) and X'^\j (|]|) can be represented in a form similar to 

(ED 



X^'p = PV-'X^\jVP (3.7) 

here V is the generator of the canonical transformation determined in the physical sub- 
space 

Vss = V2, K,. = K,a , = K,, = (3.8) 

and P is the projector on physical subspace 

P = (l/4)(2 -n)[2 -n- 2n{h+ah)s]. (3.9) 

Naturally the operators P and V commute. The partition function of the Hubbard model 
can be presented in the form similar to (p.5|), (|2.6|) 



Z = Tr(Pexp(-/5if(X^V)) = Tr{P exp{-l3H{XZi)) (3.10) 

Since the energies of the two-hole states and the triplet states are essentially higher , one 
can omit the projector P in ( 3. 10 ) at the low temperatures and work with the transformed 



Hamiltonian . This approximation is valid for small number of holes so far as the energy 
of the two-hole states at the lattice site and the triplet state are situated approximately 
at the center of the singlet hole bands. 



4 



IV. Renormalization of the hole energy in an anti- 
ferromagnetic state of the t — J model 



In this part of the work we will concern ourselves with an antiferromagnetic state for 
the t — J model [0, |, [T^ and will discuss the hole energy. 

As was recently proposed by Kane, Lee and Read (KLR hereafter) in the limit that 
the exchange energy J is much less than the hopping matrix element t, the hole spectrum 
is strongly renormalized by the interaction with spin excitations. It was established that 
the hole can be described by a narrow quasiparticle band with the quasiparticle residue 
of an order J/ 1 and a bandwidth of an order J, both for the Neel and RVB groundstates. 
Here we consider the groundstate of the local spins to be a quantum Neel state. By 
applying KLR approach we will show that our case is rather more complicated than in 
"holon" model, but the results are the same. 

An antiferromagnetic state for the Hubbard model was considered in Ref. [^. The 
Hamiltonian of the t — J model in Ref. 0] was obtained on the basis of the nonhermitian 
representation for the Hubbard operators ( pl6|) 

Ht-j = t K,[-l/2-{hlK) + {S^a)]K + J SnS„. (4.1) 

<n,n'> <n,n'> 

Retaining the notations of [0 we will write the Hamiltonian of the model. 

The Heisenberg part of the Hamiltonian ( [4.1| ) expressed through the bose-operators b 
and c, connected hj u — v transformation with the primary Dyson-Maleev bose-operators, 
which are associated with two Neel subluttices. It has the form: 

k 

a = ^J(l-7k)'/', 

7k = l/2{cos{kr,a) + cos{kya)). (4.2) 

Next we consider the addition of holes. The hopping Hamiltonian ( |4. 1|) expressed in 
terms of two-subluttice excitations and holes in lower band is 

^spl ^ 1/2^ |t^|(f/^ + V^)\{hX + C_Jp+ -qPTk + (&q + Jp^k+qPik] 
k,q 

H^^' = l/2ri' Y |tk|[pWTk(-l/2 + VlU 

k,q,k' 

t^k'f^k'+q^'k'+qV + ^k'Vk'+qCk'+qCk' + ^q-k' Vk'&lk/+qCk' + 

^-k'f^k'-qC-k'&k'-q +Pjk-qPik(^ ^ c)], (4.3) 

where r/^ = l/^z; = ((1 + u^)/{2u^)f\ = ((1 - u^)/{2u^)f\ = {1 - 7^)^^^ 
Pok'Pak are holes operators in the lower hole band. The first two terms in the round 
brackets in if gives in mean-field approximation the primary dispersion of the holes: 
= l^kl/V^, where z -number of the nearest neighbours, tk = zt'y^ ■ 
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The terms in ( [4.3|) with unconserving of the number of the particles give the two self- 
energy diagrams at zero temperature corresponding to emission and absorption one and 
two magnons. 



So the differences between the "holons", or spinless holes approach j^, |T0[ and this 
model are: (1) the existence of the primary mean-field dispersion of the hole ,(2) two 
types of the self-energy diagrams which strongly renormalized the hole spectrum. 

One can obtain from the Hamiltonian ( [4.3|) the expression of the self-energy for the 
hole with particular spin projection. 

It consists of the two contributions containing emission and absorption one and two 
magnons. Since the model is symmetric over the spin projections the self-energy in the 
non-crossing approximation (KLR) has the form: 

E(k, ^) = + ^ /i(k, q)G(k -ci,u;-E^) + 
q 

J2 /2(k, q, k')G(k -q,uj-E^,- Ek'-q), (4.4) 

q,k' 

Where /i and /2 - "two-hoppes" functions, which contain the information about dif- 
ferent types of coupling of a hole with spin excitations. One can easily get from 



/l(k,q) = (l/4)tktk-q(t/q + Vq)2 

/2(k, k', q) = Ztktk-qV^k'-qf^k'- (4.5) 

In the KLR method we have the selfconsistent integral equation for the hole propaga- 
tor: 

G(k,o;) = i— -. (4.6) 

uj — L(k, ujj 

Using the dominant pole approximation we write the hole propagator as, 

G(k, CO) = + G.„„ (4.7) 

— CJk + «I k 

where = -ReS(k, ci;k),and Fk = /mS(k, c^k)- As pointed out in Ref. [^] there is a general 
statement that Fk = for all of the low-energy poles. Since we are interested in the low 
energies, we can set Fk = 0. 

The quasiparticle residue in ( |4.7] ) 



may be estimated by omitting the incoherent part in (^4.7|) 



flk-q 



ak< [l+E/i(k,q). ^ 



E /2(k, q, kO . ^ %^ ^ (4.9) 
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By applying the scheme of we will evaluate the integrals in the denominator of 
( [4.9[ ) and will show that both integrals are of the order t/ J » 1. 

Near the bottom of the spectrum, where cufc = and in the case J = we get the 
divergence of the integrals at the small |q|: 

/i = (t^zV8v^)7^.Elql 



2m* 



q.k' 



Ok* 



2m* 



2 ' 



(4.10) 



where m* is the mass at the bottom of the band. Throughout this work the distance 
between lattice sites a = 1. The both integrals ( |4.10|) are diverged in the dimensions two 
and three. Here we will be concerned with only in the case of c? = 2. In the case of small 
but finite J we have 



flk* 



2m* 



, m Ok* 
J 



(4.11) 



Furthermore, as shown in 



m/a]^* where m depends only on t and m ~ 1/t. 



So 



t 

J' 



The integral has the more complicated form. In the case d = 2 

Ok* 

J|k'| + J|k'-q 

and its mathematical treatment is unwieldy, but the result is the same 



l2-t'J J jdH rf|k'| d^^^^^^ 



2m* 



2 ' 



(4.12) 



(4.13) 



2 m Ok* 



J 



t 

1' 



(4.14) 



thus, the quasiparticle residue Ok < J/t. 

By making physically justified assumption about the behavior of the imaginary part 
of the self-energy, one can obtain a rough estimation for Ok which is in coincides with 
( [4.14| ). Now this assumption may be checked. 

Using the Kramers-Kroning relation we have for the quasiparticle residue 

-1 



Qk = (^1 + J r(k, y)/ (y - u^fdy 



(4.15) 



As it was argued by KLR, r(k, u) vanishes like a power of {u — cuk) at {u — cjk) ^ J . At 
{uj — cjk) > J we expect that the scattering dominates and r(k, c<j) ~ t . Therefore, the 
integral in the denominator of ( [4.15| ) may be cut at y = J 

-1 



Ok 



roo 

Jj ^"^^^ 



(4.16) 
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The behavior of the imaginary part of the self-energy in non-crossing approximation 
may be obtained from expressions 

r(k,w) = -JmS(k,cj) = ri(k,cu) + r2(k,cu), 

TT 

ri(k, ^) = E /i(k, q)A(k -q,u;-E^), 
q 

r2(k, uj) = Yl /2(k, q, k')A(k - q, - - E^,_^), (4.17) 

q,k' 

where y4(k, cj) = ^IniTikjU!) is the spectral function. 

For {uj — ciJk) ^ ^, the dominant contribution to r(k, tu) will come from the pole of 
y4(k, cu) . Furthermore, the leading contribution will be given by the very small |q| and 
|k'| . For ri(k, we may write 

r,{k, u) ^ {ez'/S^2)^W E H^i-^ - o^k - J|q|) ^ I'Sf (^^)' • (4-18) 

So, ri(k, cu) ~ (cu — cuk)^ very near to the pole. Extrapolating Eqs. ( [4.17| ) in the region 
{ui — cuk) ~ J we get ri(k, Co") ~ t . For r2(k, u;) we obtain 

r2(k, c.) = (tV/2)7^^ E , ^(^ - Ik' - q| - Ik'l). (4.19) 

q,k' III ^1 

After some mathematical treatment one can get 

r2(k,a7)^— ^ , (4.20) 



J \ J 

and at (cu — uj\^ J ^ r(k, c<j) ~ t . 

Thus, we have demonstrated the application of the approach of KLR method to this 
model and have obtained the similar results, i.e. strong renormalization of a hole spec- 
trum, narrow quasiparticle band with the quasiparticle residue J/t and the bandwidth J 



V. The functional integral for the strongly correlated 
system 

Several types of functional integrals can be generated by the representation ( p.4|) for 
the Hubbard operators. We shall consider two types of such representation. The first 
can be named "constraint" representation. The second can be named as "compensation" 
representation. 
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For the first type of a functional integral we shall use some representation for spin— | 
operators 



S = 1/2(6+ (t6), As = (b+b) -1=0 (5.1) 

here b^,ba for a =],i are creation and annihilation operators which can be boson or 
fermion operators, Ag is the spin antiprojector operator which is equal to zero on the 
physical subspace and it is equal to n = —1,1,2,3 on the unphysical subspace. The 
condition = is usually named a constraint, but it is very special type of a constraint. 



The antiprojector operator Ah can also be introduce for the representation (|3.4| ) for the 
Hubbard operators: 

Ah = (l/4)(3n'^ - 2d^ + 2Sh+ah). (5.2) 

This antiprojector is equal to unity on the two-hole state and on the one-hole triplet 
states. It is equal to zero on the hole vacuum state and on the one-hole singlet state. 
The introduction of the antiprojector A is motivated by an identity connecting it with 
the projector P: 

P = (l/27r) r^exp(aA)rfA. (5.3) 
Jo 

The identity ( |5.3| ) gives us possibility to construct the functional integral for any dynam- 
ical system which can be expressed in terms of canonical variables and act on a limiting 
part of the Gilbert space which is defined by the condition A = 0. In our case of a strongly 
correlated system we have the following representation for the partition function or the 
generating functional of the Green functions 

Z{yf) = J J exp{S + I jyrJ2{Xsimk) - 1) + 



dr J \ dr 
n Dh+{r)DhUr)Db+^{r)DbMDXsi{r)DXHi{r 



S = I dTY,{hthi + btbi-H{hihuSi{btM)) + 
Jo I 

xt\hi,hMbiM))yf) (5.4) 



Here S is the action of a system; Hlh'^ , h,S{b^ ,b)) is the Hamiltonian of a strongly 
correlated system; 1 is the lattice site index; A^, Xh are the Lagrange multipliers which do 
not depend on the temperature time r what is ensured due to the presence of 6- functions. 
Because the spin subsystem possess the gauge invariance the (5(9As(r)/9r) fixes the gauge. 

Another possibility to construct the functional integral is connected with the idea of 
compensation of an unphysical contributions suggested in the work by Porov and Fedotov 
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0]. They use the equahty to zero of the Hamiltonian on the unphysical states. Let us 
add to the Hamihonian some addition 5H which is equal to zero on the physical states 
and possess the properties 

Tr{e^^{-I35H))HP = 0, [H, 6H] = (5.5) 



Summation in (5^) is produced over the unphysical states only . Then the total Hamil- 
tonian Ht = H + 6H can be used in the total Gilbert space without any constraints for 
computation of the partition function. The functional in this case has simple form: 

Z{ya,) = J exp{Sr) n Dhl{T)DhUr)Dbl{r)Dk^{T), (5.6) 

l,T,a 

here b^^, bia are the Fermi fields and the addition to the action comparing with (|5.4|) has 
form: 

6S = ~f36H = 1/2 ^[?7r(6+6; - 1) - {in + \n{3/2)){3h+hi + 2Sihlahi)] (5.7) 
I 

This addition to the action is nonhermitian. In fact we can easily construct a two- 
parametric form of 6S. 



VI. The representation of the total Hubbard algebra 



In some physical situations for example when the metal-insulator transition is studied 
it can be necessary to describe holes in the lower Hubbard band and electrons in the upper 
Hubbard band simultaneously. In this case it will be convenient to have representation 
that describes these electrons and holes as the different degrees of freedom simultaneously. 
For realization of this possibility we shall consider the additional Gilbert space consisting 
of the states of spin— | and two types of fermion with spin— | for describing electrons 
in the upper band and holes in the lower bonds. The mapping of the physical electrons 
states |0>,|t>)i>)|2> into states of the additional Gilbert space is following: 

I T> ^ \s T), I i> ^ \s i) 

|0> ^ {l/V2)[h+\s i)~hl\s^)], 

|2>^ (l/v^)[e+|si)-e+|sT)]. (6.1) 

Here \sa) for a =t, i are the states of the spin— |; /i^, /Iq,, e^, are the creation and 
annihilation the Fermi operators of holes in the lower band and electrons in the upper 
band. The physical states are singlets formed from the local spin— | and the hole (the 
electron) spin— |. One can get a representation of the total Hubbard algebra based on 
the described mapping. The representation of the lower subalgebra has form: 

X"^^ = {l/A){n^ -2d^ -2t^)Pe. (6.2) 
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The representation of the upper subalgebra is quite similar 

= 2aV2{e+S)-ail - rf)Ph 

X'''' ^{1/A){rf -2(f -2f)Ph. (6.3) 

Here we use the following short notation: 

,5= (1/4)(1-2S(7), P^^l-rf + d", 
n^ = (c+c), d" ^c^c^clci, 

t" = Sc+(7C, far c:^h,e . (6.4) 
The representation of the middle subalgebra has form: 

N ^X^^ + X^^ = P^Ph 

S = {l/2)aa/3X^^ = sPePh, (6.5) 

here s is operator of the local spin. Using this representation for the Hubbard operators 
we can easily construct the functional integral and the diagram technique. Such type 
representation can be also useful in variational approach for the construction of the explicit 
form of interaction the charge carriers with the local spins. 
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